UNCLASSIFIED 


AD  NUMBER 


ADB280586 


NEW  LIMITATION  CHANGE 
TO 

Approved  for  public  release^  distribution 
unlimited 


EROM 

Distribution:  Further  dissemination  only 
as  directed  by  Office  of  Scientific 
Research  and  Development^  Washington^  DC 
20301,  DEC  1941,  or  higher  DoD  authority. 


AUTHORITY 


OTS  index  dtd  Jun  1947 


THIS  PAGE  IS  UNCLASSIEIED 


rvt>f?c 

f^-AjL 

Lu  . 


Copy  No. 


OS/t  0  6 


INTERNAL  BALLISTICS  OF  POWDER  DRIVEN  ROCKETS 


DISTRIBUTION  STATEMENT  F: 

Further  dissemination  only  as  directed  by 

_ Cp/re  ii4 

or  higher  DoD  authority. 


fey 


OSRD 

Washington  DC  20330 


Emory  Lakatos 


>ry  .-v  J  -f-'n'"'*-  n  —  •'  « 


Reproduced  From 
Best  Available  Copy 


( 

?* 


20020701  205 


NATIONAL  DEFENSE  RESEARCH  COMIITTEE 


REPORT  NO.  A-22  :  PROGRESS  REPORT 


INTERNAL  BALLISTICS  OF  POiVLER  DRIVEN  ROCKETS 


by 

Emory  Lakatos 


Approved  December  16,  19Ul 

for  submission  to  the  Section  Chairman 


Approved  December  18,  I9I4.I 

for  submission  to  the  Division  Chairman 

C.  N,  Hickman, 
Chairman,  Section  H, 
Division  A 

Approved  December  I8,  I9I4I 

for  submission  to  the  Committee 


^ o  Cil  qA^  (  T^u- 1 
Richard  C,  Tolman 
Chairman,  Division  A 


Emory  i'Lakaios,  Author 
Consultant,  Section  H, 
Division  A 


Preface 


This  report  is  pertinent  to  the  project 
designated  by  the  War  Department  Liaison 
Officer  as  OD-26.  It  pertains  to  work  in 
progress  at  Indian  Head,  Maryland,  as  a  joint 
project  of  the  Bureau  of  Ordnance,  U.S.  Navy, 
and  Section  H,  Division  A,  National  Research 
Committee. 


Distribution  of  copies  of  thio  report.  — 

The  report  was  prepared  for  duplication  bn 
December  20,  19U1.  The  initial  distribution 
of  copies  was  as  follows; 

Copies  No,  1  to  23,  inclusive,  to  the  Office 
of  the  Secretary  of  the  Committee  for  distribu¬ 
tion  in  the  usual  manner; 

Copies  No,  2k  and  25  to  the  Chief  of  the 
Bureau  of  Ordnance,  Navy  Department,  Washington, 
D,C,  (Attention;  Research  and  Development  Divi¬ 
sion)  , 

Copies  No,  26  and  27  to  the  U,S,  Naval  Powder 
Factory,  Indian  Head,  Maryland  (Attention: 
Lieutenant  Commander  H,  J,  Orth) ; 

Copy  No,  28  to  Major  L,  A,  Skinner,  Ordnance 
Department,  U,  S,  Array; 

Copy  No,  29  to  C,  N,  Hickman,  Chairman, 
Section  H,  Division  A; 

Copy  No.  30  to  Dr,  Vif,  A.  Fowler,  Consultant, 
Section  H,  Division  A. 

Copy  No,  31  to  Prof,  E,  C,  Watson,  Member, 
Section  H,  Division  A. 

Copy  No.  32  to  E.  Lakatos,  Consultant, 

Section  H,  Division  A. 


CONTENTS 


Section  Page 

Abstract  . 1 

1.  Basic  equation  . . . 2 

2 .  Law  of  burning . 2 

3.  Additional  relationships  .  U 

U.  Equilibrium  pressures  .  5 

Variation  of  pressure  virith  S/A;j^,  the 
ratio  of  burning  area  to  throat  area  ...  8 

6.  Variation  of  pressure  with  fraction  of 

propellant  consumed  .  9 

7.  Exit  velocities  of  gas  in  simple  con¬ 

vergent  nozzles  and  in  expanding 

nozzles  . 12 

8.  Exit  area  .  17 

9.  Comparison  of  thrusts  obtained  with  a 

simple  convergent  nozzle  and  with  an 

expanding  nozzle  for  the  case  of 

constant  pressure  .  18 

10.  Projectile  velocity  .  19" 

1 1 .  Efficiency  for  the  case  of  a  constant 

chamber  pressure  . . . 22 

12.  Muzzle  velocity  expressed  in  terms  of 

the  properties  of  the  propellant  .  2$ 

13*  Rocket  design  calculations  .  25 

lU.  Performance  of  geometrically  similar 

shells  .  37 

Appendix .  Ul 


V 


List  of  Figures 


Figure  Page 


1.  Effect  of  finite  pressure  ratios  on  the 

gas  exit  velocity  from  an  expanding 

nozzle . .  13 

2.  Ratio  of  nozzle  mouth  area  to  nozzle 

throat  area  [Eq.  (22a)]  .  lli 

3.  versus  (Vq/V-)Z  for  several  pairs 

of  values  of  and  A  /X  . . .  l5 

U.  Ratio  of  burning  area  to  throat  area 
versus  equilibrium  chamber  pressure, 
for  two  different  lots  of  powder  .  I6 

5.  Propulsion  efficiency  versus 

X't^Mym^]  .  29 

6.  Muzzle  velocity  as  a  function  of  pro¬ 

pellant  charge  and  pressure . 30 

7.  TNT  load  versus  propellant  charge  and 

pressure  . .  30 

8.  Effective  jet  velocity  and  thrust 

coefficient  versus  pressure,  for  two 
different  weights  of  propellant  charge  . .  31 

9*  Burning  time  versus  pressure  expansion 

ratio  . . . .  31 

10.  Shell  structure  . .  32 


vi 


INTERNAL  BALLISTICS  OF  POVffiER  DRIVEN  ROCKETS-'^ 


Abstract 


An  approximate  theory  of  the  internal  ballistics  of  rockets 
driven  by  colloidal  propellants  is  developed  in  this  report.  A 
number  of  approximate  relationships  are  derived  which  may  be  of 
use  in  preliminary  design  work.  Particular  emphasis  is  placed 
on  the  case  of  a  constant  chamber  pressure.  Included  in  the 
report  are  sample  design  calculations  for  a  3“in.  l5-lb  shell. 
The  calculations  are  summarized  on  graphs  which  show  the  weight 
of  TNT  that  can  be  carried  as  a  function  of  the  propellant 
charge  weight,  chamber  pressure  and  muzzle  velocity.  The 
results  show  that  the  largest  weight  of  TNT  can  be  carried  when 
the  operating  pressure  is  between  25  and  50  atmos.  This  result 
is  then  shown  to  apply  equally  well  to  all  shells  geometrically 
similar  to  the  one  considered. 

The  theory  developed  in  this  report  is  necessarily  approx¬ 
imate  in  the  following  respects; 

(a)  In  formulating  the  law  of  burning,  heat  losses  in 
walls  are  neglected.  This  means  that  no  account  is  taken  of 
the  fact  that  the  mean  temperature  of  the  gas  in  the  rocket 
chamber  is  less  than  the  temperature  of  the  layer  of  gas  im¬ 
mediately  adjacent  to  the  burning  surface. 

(b)  The  mass  of  gas  discharged  through  the  nozzle  of  the 
rocket  is  calculated  on  the  assumption  that  the  steady  flow 
formulas  may  be  applied  to  variable  flow  problems  on  a  quasi¬ 
steady  state  basis.  A  more  rigorous  treatment  of  this  problem 
does  not  seem  to  be  available  at  this  time. 


1_/  This  report  is  a  revision  of  one  of  October  29,  19l;0, 
For  other  material  pertinent  to  the  present  subject,  see  C.  N. 
Hickman,  MDRC  Report  A-k,  Appendix  D;  C,  N.  Hickman, 

Memos.  A-18M  to  A-22M;  J.W.M.  Dulfond,  NDRC  Report  A-21|.  The 
last-named  report  is  a  relatively  elementary  treatment  of  the 
mechanical  efficiency  of  rockets. 
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(c)  The  variation  in  gas  temperature  in  the  initial  stages 
of  burning  is  ignored,  Ihile  this  simplification  yields  a  cor¬ 
rect  prediction  of  the  equilibrium  pressure,  it  may  at  times 
lead  to  incorrect  results  in  the  prediction  of  how  soon  equilib¬ 
rium  is  reached. 

(d)  Effects  of  friction  are  ignored. 

The  English  system  of  units  based  on  the  pound  as  the  unit 
of  force  and  slug  as  the  unit  of  mass  is  used  throughout  the 
report.  The  list  of  symbols  used  is  given  in  the  Appendix. 

T.  Basic  equation 

The  basic  equation  connecting  the  rate  of  burning  and  the 
rate  of  gas  discharge  is  that  the  mass  burned  up  to  time  t  is 
equal  to  the  mass  of  gas  in  the  chamber  plus  the  total  mass  of 
gas  discharged.  An' equivalent  statement  is  that  the  rate  at 
which  the  propellant  is  consumed  equals  the  rate  at  which  the 
gas  is  discharged  plus  the  rate  at  which  the  mass  of  gas  inside 
the  chamber  changes j  in  symbols, 

m  »  m^  +  (dm^dt),  (1) 

where'®  is  the  mass  (slugs)  of  propellant  burned  per  second, 
fflj  is  the  mass  (slugs)  of  gas  discharged  per  second  and  m  is 
the  mass  (slugs)  of  gas  in  the  chamber  at  the  time  t. 

2,  Law  of  burning 

From  the  outset  vire  shall  restrict  the  treatment  to  the  case 
of  constant  burning  surface,  such  as  a  tube  burning  on  its 
inside  and  outside  curved  surfaces,  or  a  solid  cylinder  burning 
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only  on  its  flat  ends.  Let  D  be  the  dimension  of  the  propellant 
along  the  line  of  burning.  This  quantity  is  so  defined  that  I'D 
is  the  greatest  depth  burned  away  below  any  receding  surface 
before  the  whole  grain  is  consumed.  In  the  case  of  a  tube,  D 
is  the  wall  thickness.  For  a  cylinder  burning  from  both  of 
its  flat  ends,  D  is  the  length  of  the  cylinder.  If  the  cylinder 
burns  from  only  one  end,  D  is  taken  as  twice  the  length  of  the 
cylinder.  It  is  easy  to  see  that,  in  the  cases  cited,  the 
definition  is  satisfied.  For  example,  in  the  tube,  the  wall 
burns  a  distance  ^  from  the  inside  and  a  distance  gD  from 
the  outside. 

With  this  understanding,  we  take  as  the  law  of  burning, 

B  ^  =  b  +  BP,  (2) 

where  Z  is  the  fraction  of  the  propellant  consumed  at  time  t, 
b  and  B  are  constants  of  the  specific  propellant  used,  and  P 
is  the  chamber  pressure  (Ib/ft^).  The  physical  meaning  of 
Eq,  (2)  can  be  shown  as  follov/s. 

If  it  is  assumed  that  the  burning  takes  place  by  parallel 
layers  and  that  the  total  burning  surface  is  constant,  the 
fractional  mass  of  propellant  consiimed  will  be  proportional  to 
D*Z,  the  total  depth  burned  away.  Hence  D(dZ/dt)  must  be  the 
total  distance  that  the  flame  travels  per  second.  That  is  to 
say,  if  the  flame  velocity  is,  say,  1  ft/sec,  then  D(dZ/dt) 
will  be  2  ft/sec  for  the  two  burning  surfaces.  It  is  also 
clear  that  b  has  the  dimensions  of  a  velocity  (ft/sec)  and 
that  B  is  a  velocity  per  unit  of  pressure  [ (ft/sec)/(lb/ft^) ] . 


-k- 


It  should  be  remarked  that  this  form  of  the  burning  law 

has  been  strongly  criticized  in  the  past  because  it  is  not 

2/ 

correct  either  physically  or  philosophically.-'  However, 
the  writer  has  found  that  it  appears  to  be  the  only  simple 
relationship  using  the  burning  constants  detemined  from 
closed-bomb  experiments  that  correctly  predicts  the  equilib¬ 
rium  pressure  in  rockets.  For  example,  the  theoretical  law 

2/ 

derived  by  Crow  and  Grimshaw-'  fits  the  experiments  with  bombs 
but  for  rockets  gives  materially  Incorrect  values  of  the 
equilibrium  pressure  attained  as  a  function  of  the  ratio  of 
burning  area  to  throat  area.  Just  why  apparently  equivalent 
laws  should  give  different  results  is  not  yet  clear  to  the 
writer. 

3.  Additional  relationships 

The  remaining  equations  necessary  for  the  formulation  of 
the  problem  are 

m  =  m^dZ/dt,  (3) 

where  is  the  initial  mass  of  the  propellant.  The  gas 

equation,  assuming  that  the  gas  can  be  regarded  as  an  ideal 
gas,  is 

Pv  =  aj,.  (U) 

where  P  is  the  pressure  (Ib/ft*),  v  is  the  specific  volume 
(ft^/slug),  R  is  the  gas  constant  (R  =  1930  ft  Ib/slug  °F) 
and  T  is  the  combustion  temperature  {'^F  absolute). 

''  “  '  .f"  '  ■  ^ 

2/  In  particular,  see  Crow  and  Grimshaw,  Phil.  Trans, 

Roy, -Soc.  (London)  23Qj  A691 ,  pages  387-Un  (1S\32). 
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Also,  the  instantaneous  value  of  the  mass  of  gas  in  the 
chamber  is  given  by 

mg  =  VA,  (5) 

where  V  (ft^),  the  instantaneous  gas  volume  in  the  chamber, 
may  be  expressed  as 

V  =  Vq  +  (m^/,9  Z  =  +  Vp-  Z,  (6) 

Vq  being  the  initial,  or  clearance,  volume  (ft^)  V  ^the  volume 
of  the  propellant  (ft^)  and  6  the  density  of  the  propellant 
(slug/ft^). 

Finally,  the  rate  of  gas  discharge  may  be  written  as 

m^  =  kP,  (7) 

if  the  chamber  pressure  does  not  fall  below  about  ]i300  Ib/ft^, 
or  2  atraos,  above  zero  pressure.  The  constant  k  has  the  value 


,r? 


“  A/rt 


2 

+  1  - 


(«  +  -  1) 


(8) 


or 

k  =  (h  A^,  (8a) 

where  is  the  throat  area  of  the  nozzle  (ft®),  X  is  the 
adiabatic  constant  —  that  is,  the  ratio  of  the  specific  heats 
at  constant  pressure  and  volume  —  and  the  factor  ,  which 
represents  the  quantities  under  the  radical  in  Eq.  (8),  is  a 
constant  that  depends  on  the  properties  of  the  specific  pro¬ 
pellant  used. 


U.  Equilibrium  pressures 

By  combining  Eqs.  (U),  (5)  and  (6)  we  can  write  the  expres¬ 
sion  for  the  mass  of  gas  in  the  chamber  as 


(9) 
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Equation  (1)  may  then  be  written  as 

®  RT  dt 

*  Id  d2 

~  ™d  ■*' RT  dZ  dt 


or,  upon  eliminating  dZ/dt  by  means  of  Eq.  {3>)  > 

m  =  Ki.  +  (1/RT)  V  +  (V  +V  Z)  (m/m,). 

P  .0  p  dZ  J  ^ 

Solving  for  the  derivative  dP/dZ,  we  obtain 

RltaJI  -  (m^/m)]  -  P  V 


dZ 


^o  ^  \  2 


(10) 


For  equilibrium,  dP/dZ  must  be  zero.  The  equilibrium  pressure 
is  therefore  given  by 

RTm^  mo 


But  ni,/V_  is  the  density  ^  of  the  propellant,  and  RT^  is 
X  p  - 

the  bomb  pressure  that  would  exist  with  100  percent  density 

of  loading  if  the  temperature  were  that  for  burning  at  constant 

pressure.  Let  this  pressure  be  denoted  by  P  .  Then 

b 

P3  -  Fj,  (1  -^).  (11) 

-  m 

This  result  is  importah't  and  requires  some  comment.  First 
of  all,  it  holds  independently  of  any  particular  law  of  burning. 
It  is  true  whether  the  temperature  transient  is  taken  into 
account  or  not.  It  shows  that,  for  an  equilibrium  pressure  to 
exist,  the  rate  of  discharge  must  be  less  than  the  rate  of 
burning,-  Finally,  it  shows  that  the  value  of  the  equilibrium 
pressure  is  apt  to  be  critical  as  to  factors  that  affect  either 


1 
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the  burning  or  the  discharge  rates.  This  is  so  because  is 

of  the  order  of  200,000  Ib/in.^,  and  the  value  of  P  desired 

s 

for  rocket  work  hardly  ever  reaches  10,000  Ib/in, ^ 

If  the  burning  law,  Eq.  (2),  and  discharge  rate  law, 

Eq.  (7),  are  substituted  into  Eq.  (10),  there  results 


1 


k  P. 


m  B 

^  (|  *  t>e) 


=  0. 


Let 


and 


m^S/D  =  C, 

(12) 

b/BP^  =«, 

(13) 

The  equation  that  gives  the  (dimensionless)  equilibrium  pres¬ 
sure  ratio  explicitly  is 

A  2  ^  X  [  (k/C)  +<:*>-  1  ]  -  =  0.  (lU) 

Instead  of  solving  for  A^,  if  is  more  convenient  to  solve  for 
k/C  in  terms  of  the  equilibrium  value;  thus 


k/C  =  (o^/z^lg)  +  1  -  OC  -  Ag.  (15) 


If  S  is  the  total  burning  area, 

m  =  i  S  SD* 
z 


Then 


k/C  =  2/'iA^/^SB, 


(l5a) 


or 


\  _  ilf 

S  "  2/J 


(16) 
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5.  .Variation  of  pres Sure -with  ratio  of  burning, area 

to  throat  area 

To  see  how  these  considerations  work...  out,  consider  a 
propellant  particularly  well  adapted  for  rocket  work;  namely, 
one  having  a  composition  of  UO  percent  nitroglycerine  and 
60  percent  nitrocellulose,  without  any  stabilizer,  and  having 
properties  approximately  as  follows; 
i  =  1.25, 

R  =  1930  lb  ft/slug  °F, 

T  =  5050°F  absolute,  at  constant  pressure, 

b  =  0.065  ft/sec, 

B  =  U.5  X,  lo”"^  (ft/sec )/(lb/ft^), 

^  *  3.10  slug/ft^  [=  100  lb  of  mass/ft^]. 

Then,  from  Eq.  (U),  *  3.02  x  10 Ib/ft^,  or  210,000  Ib/in.®; 

from  Eqs.  (8)  and  8(a),  fh  -  11.97  x  10  ^  slug/lb  see,;  from 

-3 

Eq,  (13),  cx.  =  h.80  X  10  . 

By  assuming  successive  values  of  that  is,  of  Pg/P^^j 

and  substituting  in  Eq,  (16),  we  get  the  curve  in  Fig.  U 
showing  the  equilibrium  chamber  pressure  in  terms  of  the  ratio 
s/a  .  The  circles  and  squares  show  experimental  points  for 

U 

two  batches  of  powder.  The  agreement  with  the  theoretical 
curve  is  quite  good.  It  will  be  noted  that,  at  low  pressures, 
the  rate  of  change  of  pressure  with  area  ratio  S/A,j.  is  small. 
However,  at  the  higher  pressures,  the  actual  pressure  is  quite 
sensitive  to  comparatively  small  changes  in  the  area  ratio. 
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f 


t 


6.  Variation  of  pressure  with  fraction  of  propellant  consumed 
Returning  now  to  the  general  equation  [Eq.  (10)]  we  can 
rewrite  it  in  the  dimensionless  form 


(VoAp)  *  2 


X»  ■*  A[(ic/c)  *  -  D- 


This  rather  formidable  looking  equation  purports  to  show  how 
the  actual  pressure  approaches  the  equilibrium  pressure  as  a 
function  of  the  fraction  burned.  It  will  be  noticed  that  if 
—  that  is,  if  the  black  powder  charge  starts  the 
combustion  off  at  just  the  right  pressure  —  the  pressure 
throughout  the  burning  interval  will  be  the  equilibrium  pres¬ 
sure.  For  other  values  of  the  initial  pressure,  the  results 
are  more  complex.  Figure  3  shows  Ag  plotted  as  a  function 
of  (Vp/V^)Z  for  pairs  of  values  of  and 

It  will  be  noted  that  the  curves  of  Fig.  3  have  a  dual 
significance.  The  abscissa  is  (Vp/VQ)Z.  Hence  for  Z  =  1 
the  curves  give  the  pressures  attained  for  various  densities 
of  loading  at  the  end  of  the  combustion  period.  For  Z  <  1  and 
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cpnstant  loading  density  the  curves  show  the  manner  in  which 
the  pressure  builds  with  the  fractional  mass  of  propellant 
consumed.  For  example,  consider  the  curve  marked 
Xg  »  0,01,^ .  Xq/ This  means  that  we  are  considering 
the  case  where  the  equilibritun  pressure  is  0,01  x  210,000,  or 
2100  Ib/in,^  Suppose  the  design  is  such  that  Vp/V^  =2, 
which  corresponds  to  a  density  of  loading  of  66,7  percent. 

Then,  since  the  propellant  is  all  consumed  when  Z  =  1 ,  we  see 
that,  for  (Vp/Vq)Z  =  2,  X/Ag  *  0.7,  or  the  pressure  reached 
is  114.70  lb/in.®j  that  is,  the  equilibrium  pressure  is  not 
reached  (according  to  this  equation)  during  the  combustion 
process.  If  we  wish  to  know  the  pressure  when  the  combustion 
is  30  percent  complete,  Z  *  0.3,  and  (Vp/VQ)Z  =  0.6.  The 
ordinate  is  0.605  ,and  hence  the  pressure  reached  is 
0.605  X  210,000,  or  1270  Ib/in.^  If  the  density  of  loading 
is  increased,  the  same  curve  shows  that  the  final  pressure 
approaches  the  equilibrium  value  of  2100  Ib/in.^,  and  for 
100  percent  density  of  loading  —  that  is,  Vp/v^  =*  co  — 
equilibrium  is  reached  at  the  very  end  of  the  process. 

For  lower  values  of  the  equilibrium  pressure  >  say 
Xg  *  0,0025,  or  a  pressure  of  525  Ib/in,^,  the  approach  to 
equilibrium  tends  to  be  more  rapid j  but  it  is  still  not 
completed.  A  similar  situation  is  predicted  if  the  initial 
pressure  is  higher  than  the  equilibrium  value. 

These  predictions  are  only  partially  confirmed  by  experi¬ 
ment.  .?/hile  the  evidence  is  still  not  complete,  it  appears 
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that  in  most  practical  cases,  the  pressure  builds  up  to  the 
equilibrium  value  quite  rapidly  and  stays  there  until  the  com¬ 
bustion  is  completed.  By  a  "practical  case"  is  meant  one 
where  the  density  of  loading  exceeds,  say,  2$  percent  and  the 
burning  time  exceeds,  say,  O.O3  sec.  Under  these  conditions 
the  equilibrium  pressure  is  quickly  reached  even  for  very  low 
starting  pressures. 

On  the  other  hand,  if  the  burning  is  very  rapid,  lasting, 
say,  0.01  sec,  the  equilibrium  pressure  is  never  reached.  Also, 
if  the  density  of  loading  is  small,  equilibrium  is  not  attained 
even  when  the  combustion  lasts  as  long  as  O.3  sec.  The  experi¬ 
mental  curves  appear  quite  similar  to  those  in  Fig. .3,  even  in 
respect  to  the  effect  pf  initial  pressure. 

In  view  of  the  only  partial  agreement  between  experiment 
and  the  simple  theory  used  here  and  the  lack  of  a  criterion  of 
applicability,  we  are  forced  to  conclude  that  the  theory  cannot 
be  relied  on  to  predict  the  transient  state.  Hovrever,  it  does 
appear  to  predict  equilibrium  values  with  a  reasonable  degree 
of  accuracy.  Hence,  for  design  work,  a  simple  theory  of  this 
type  is  adequate  fop  most  purposes  when  tempered  by  a  knowledge 
of  experimental  results. 

It  may  be  remarked  here  that,  for  the  case  of  constant 
pressure,  the  burning  time  T  may  be  calculated  from  Eq.(2). 

This  equation  becomes,  for  this  case. 


T  =  D/(b  +  BP) 


(18a) 


12  - 


In  design  calculations,  once  the  throat  area  is  knovvn,  the 
time  T  may  be  algo  calculated  from  the  equation 

T  »  m^/m^  =  m^//i  A^P. 

7.  Exit  velocities  of  gas  in  simple  convergent  nozzles  and 


3/  J.  Franklin  Institute  230,  U33~U5l  (I9ii0). 


AOi 


A  fee/"  of  pressure  rofos  on  f/ye  gas  ex  A?"  \re/oc/fg  /^ro/77  an 

expan c//nq  nozz/e. 


/7^  3.  A/ Ae  /cr severe/ pa/rs  of 

of  Xe  and  Xo^  X^. 


J 


4.  /f?7//c7  of  bi/rn/ncf  oreo  bo  bdroo/  area 
yerjas  efo/Z/br/o/r)  c/?(7/7?6&r  prejsore^  for  Ztft> 
cZ/f f erenb  /oAs  of  powc/er(oya\ 
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For  the  propellant  already  referred  to  in  Sec.  5}  these 
gas  velocities  are 

=  3280  ft/sec  (19a) 

and 

Note  that,  viith  this  propellant,  the  highest  velocity  possible, 
(with  infinite  pressnre  ratio)  is  9890  ft/sec.  The  effect  of 
finite  pressure  ratios  can  be  seen  from  Fig.  1 . 

It  should  be  also  remembered  that  the  mass  of  gas  dis¬ 
charged  per  second  is  given  by  Eq. (?)  lor  both  types  of  nozzles. 
Moreover,  the  pressure  at  the  throat  for  both  types  is  given  by 

Pt  *  PI-FTt)  .  (21) 

or 

=  0.995  P  (21a) 

for  the  propellant  here  considered. 


8.  Exit  area 

The  standard  formula  for  the  nozzle  mouth  area  A  is- 

m 


which,  for  the  specific  propellant  here  considered,  becomes 


h/  For  the  theory  of  expanding  nozzles  see,  for  example, 
Martin,  Textbook  of  mechanics,  vol.V,  ''Thermodynamics",  pp, 230-236. 
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The  plot  of  this  function  appears  in  Fig^  2.  One  should  note 
that,  in  actual  designs,  i^i  is  frequently  impracticable  to  use 
full  expansion.  That  is,  is  really  fixed  by  the  diameter 
of  the  projectile  rather  than  by  theoretical  considerations. 

In  such  cases.  Fig.  2  can  be  used  to  find  the  pressure  drop 
when  the  areas  are  prescribed. 


9,  Comparison  of  thrusts -obtained  with  a  simple  convergent 
nozzle  and  with  an  expanding  nozzle  for  the  case  of 
constant  pressure 

Assume  that  the  orifice  area  of  the  simple  nozzle  is  the 
same  as  the  throat  area  A^  of  the  expanding  nozzle.  The 
thrust  F^  (lb)  is  the  time  rate  of  change  of  momentum 
(slug  ft/sec2)  plus  the  imbalanced  force  (lb)  acting  on  the 
orifice  area,  or 

f't  ■  *  h  <23) 

For  practical  purposes  the  atmospheric  pressure  P.  may  be 

Or 

neglected  in  comparison  with  the  pressure  P^,  If  we  then 
substitute  the  value  of  from  Eq,  (7),  w^  from  Eq.(19),  and 
P^  from  Eq.(21),  we  get 


2  \l/(^-1) 


+1 


(2i4) 


For  the  propellant  cited  (Sec,  5), 

»  1.2U  A^_P,, 


(2iia) 


The  numerical  factor  in  Eq,  (2lia)  is  called  the  thrust 


coefficient.  It  is  a  very  convenient  criterion  for  making 
comparisons  in  performance,. 
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In  the  case  of  an  expanding  nozzle  of  correct  ratio,  that 
is,  one  which  expands  the  gas  to  atmospheric  pressure,  the 
thrust  is  given  by  the  first  term  of  Eq. (23),  except^  of  course, 
that  the  value  of  the  exit  velocity  is  now  w^^^,  not  w^.  Suo- 
stituting  the  value  of  w^^  from  Eq.  (11) ,  we  get 


Fg  -  PS  j  3  +1  j 

Numerically, 


2  ^(Jf +!)/(« -1)  /  "PfJIfTjTT 


j  ’  "I Pi 


.  (25) 


r 


2.03  A.P 


y.-ipi 


0.200 


(25a) 


By  referring  to  Fig.  1  we  find  that,  for  a  pressure  of 
200  atmos,  the  radical  in  Eq.(25a)  is  0.81.  In  this  particu¬ 
lar  case,  the  thrust  coefficient  would  be  0.8l  x  2,08,  or  1*69* 
Comparing  this  value  vfith  that  for  the  simple  orifice  [Eq.(2ha)] 
we  see  that  the  thrust  is  improved  in  the  ratio  of  1,69  to  1.2l|, 
or  1.36,  For  lower  chamber  pressures,  this  improvement 
decreases,  owing  to  the  decreasing  value  of  the  radical. 

For  the  case  of  an  overexpanding  nozzle,  reference  should 
be  made  to  Malina's  paper. The  method  to  be  used  in  the  case 
of  under expans ion  is  illustrated  in  Sec,  13* 


10.  Projectile  velocity 

To  avoid  the  confusion  that  sometimes  may  arise,  we  will 
base  the  momentum  calciilations  on  first  principles. 


5/  See  reference  3* 
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Consider  a  rocket  of  instantaneous  mass  Mp(t)  moving  to 

the  right  and  having,  at  time  t,  a  velocity  x  relative  to 

the  ground.  Take  i  as  positive  in  the  right-hand  direction. 

An  observer  stationed  on  the  rocket  causes  the  release  of  a 

differential  mass  of  gas  dm^  with  a  velocity  w^,  this  velocity 

being  measured  relative  to  the  rocket  after  the  mass  dm,  has 

-  d 

been  released.  The  observer  takes  the  origin  of  his  coordinate 

frame  as  being  on  the  rocket.  Relative  to  this  frame  he  calls 

positive  when  it  is  in  the'  left-hand  direction. 

Since  it  is  assumed  that  no  external  forces  are  acting, 

the  law  of  conservation  of  momentum  requires  that  the  momentum 

before  dm^  is  released  is  equal  to  the  algebraic  sum  of  the 

momentums  of  the  rocket  and  of  dm  after  the  latter  is  re- 

d 

leased.  In  other  words,  if  the  velocity  of  the  rocket  after 
the  release  of  dm^  is  x  +  dx,  the  law  of  conservation  of 
moment iim  requires  that 

Mp(t)k  *  [Mp(t)-dm^]  [x+dx]  +  [x+dx  -  w^]dm^^ 
or 

0  »  Mp(t)dx  -  (26) 

Since  all  the  quantities  in  this  equation  are,  in  general, 
functions,  of  the  time, 

X  =  w  dm,/M  (t) , 

Jo  ^  d  p 

or 

• 

X  *  I  w  m  dt/M  (t),  (27) 

Jo  m  d  p  ^  ' 

where 

M  (t)  =  M  ~  1  m  dt, 

P  o  Jo 


(28) 
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M  being,  of  course,  the  initial  mass  of  the  loaded  projectile, 

0  - 

Formally,  at  least,  the  problem  of  calculating  the 

velocity  is  solved.  Actually,  in  the  general  case,  the  function 

w  is  of  such  nature  that  any  expression  containing  it  is  not 
m 

likely  to  be  integrable  by  any  except  numerical  methods. 

However,  for  the  special  case  of  a  constant  chamber  pres- 


sure,  the  integration  is  quite  simple  because  vf  and  m  are 

m  d 

both  constants.  For  this  case^then. 


w  m ,  I 
m  d 


- ; -  =  w  log  - ;; —  .  (2?) 

M  -m-t  M  '•m  t 

o  d  0  d 


When  the  combustion  ceases,  at  t  =  T  ,  the  chamber  is  still 
filled  vjlth  gas,  which  then  escapes  at  a  constantly  diminish¬ 
ing  velocity.  If  the  momentum  of  this  residue  is  neglected, 
we  may  assume  that 

m  T  =  m^.  .  (30) 

d 

Making  this  substitution  in  Eq.  (29),  we  obtain  for  Xj^,  the 

maximum,  or  muzzle,  velocity  of  the  projectile, 

^ — •  (31) 

M  -m^ 
o  ^ 


The  follovdng  relationships  are  sometimes  useful; 


(a)  If  is  the  mass  of  the  shell  without  propellant. 


xi/w 


m_  =  M  (e  ~  1 )  • 

X  00  ' 


(b)  If  the  ratio  m^  is  sufficiently  small, 


• 

a 

+  •  •  • 

J 

(33) 


-  22  - 


For  this  condition,  ths  acceleration  of  the  projectile  is 
approximately  constant  during  the  burning  period. 

11 .  Efficiency  for  the  case  of  a  constant  chamber  pressure 

The  expression  for  the  efficiency  is 

Kinetic  energy  of  projectile  at  muzzle 

Efficiency  =  - -  , 

rieat  energy  of  propellant 


or 


E  =  |[M  -m  ]xi  /J  H 


m 


(3k) 


where  J  is  the  mechanical  equivalent  of  heat,  which  is  equal 
to  776  ft  Ib/Btu,  and  H  is  the  heat  of  combustion,  expressed 
in  British  thermal  units  (Btu)  per  slug  of  propellant.  But, 
according  to  Eq.(29) 

M 


Xi^  =  VI-  " 

^  m 


log 


Li  , 

o 


and,  according  to  Eq.(20)’, 


w 


2  _ 


m 


2  ^  R  T 
^  -  1 


/ 


1 

'  P  / 


hence 


E  =■ 


^  R  T 


J  H 


_  ip/ 


/!o. 

\  m_ 


*-  1/, 


log 


1 

(Mo/m^)-1  i  • 


But  the  gas  constant  expressed  in  heat  units  —  that  is, 
R/J  (Btu/slug  ^F)  is  merely  the  difference  of  the  mean 
specific  heats  at  constant  pu'essure  and  at  constant  volume, 
provided  the  gas  may  be  treated  as  ideal;  in  other  words, 

R/j  -  Cp  -  d^. 


(35) 


(36) 


-  23  - 


Moreover,  since  most  of  the  combustion  takes  place  under  constant 
pressure,  and  in  view  of  the  simplification  it  introduces,  the 
equation 

H  =  c  T  (37) 

P 

is  taken  to  represent  the  complete  process.  Combining 
Eqs.(36)  and  (37)}  we  have 

(38) 


RT/JH  =  (Cp,  -  c^)/-C^^  =  (K  -1)/J^. 


Let 


M  /m  =  (X ' . 
o  X 


Then  Eq.(33)  becomes 


E  = 

/  p  .  (^-1)A  1 

1  - 1  - 

(»'-l) 

log  ^ - 

i 

L  1  P  / 

cc  '-1  _ 

(39) 


The  form  of  the  functions  P  /P  and  oL ’  are  shown  in 

m 


Figs.  1  and  5,  respectively.  The  latter  curve  shows  that  the 
best  value  of  ot '  is  1.25,  v^hich  means  that,  of  the  projectile 
weight,  80  percent  is  propellant  and  20  percent  is  shell. 

Assuming  an  infinite  expansion  ratio,  we  see  that  the  greatest 
efficiency  possible  is  0.65. 

As  a  more  practical  example,  consider  a  l5-lb  projectile, 
with  2.5  lb  of  propellant  charge  operating  at  200  atmos  pressure. 
Thenoc.  -  5,  X  -  1.25,  1  -!  0.65  and. 


from  Eq.(39)}  E  =  0.65  ^  0.168  =  O.lO^j  that  is,  the  effi¬ 
ciency  is  almost  11  percent. 


It  should  be  noted  that  Eq.(39),  the  expression  for  the 

efficiency,  can  be  split  into  two  factors.  One  depends  only 

on  the  pressure  ratio  and  is  the  thermal  efficiency)  it 

measures  th?  fraction  of  the  heat  energy  that  goes  into  giving 

kinetic  energy  to  tke  jet.  The  second  factor^ which  depends  only 

on  oc’,  the  ratio  of  the  initial  total  weight  of  the  projectile 

to  the  weight  of  the  propellant,  is  the  propulsion  efficiency) 

it  measures  the  fraction  of  the  jet  energy  that  is  converted 

into  kinetic  energy  of  the  rocket  itself. 

Equation  (39)  is  quite  interesting  in  that  it  shows  that 

when  the  weights  of  the  propellant  and  of  the  projectile  arfei 

fixed,  nothing  can  really  be  done  about  increasing  the  effi- 

6/ 

ciency.— '  The  value  of  X  varies  but  slightly  for  various 

proportions  of  combustion  products.  The  same  thing  is  true  for 

C  ^  )  /  X 

the  factor  1  -  (Pjjj/P)  >  foi*  the  pressure  ratios  in  common 

use.  EJven  H,  the  heat  of  combustion  of  the  propellant,  drops 
out  of  the  efficiency  expression.  This  is  not  to  say  that  the 
heat  of  combustion  H  is  unimportant...  As  will  be  shovm  in  the 
next  section,  the  muzzle  velocity  per  unit  mass  of  propellant 
varies  as  the  square  root  of  H,.  Hence,  all  other  factors 
being. equq^,  the  Substitution  of  propellants  having  larger 
heats  of  combustion  will  result  in  higher  values  of  the  muzzle 

6/  It  is  true  the  propulsion  efficiency  can  be  increased 
by  giving  the  rocket  an  initial  velocity.  Whether  the  overall 
efficiency  (including  the  energy  required  to  Impart  the  initial 
velocity)  is  also  increased  is  rather  problematical. 
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•  2 

velocity.  But,  since  varies  as  H,  no  higher  efficiencies 
can  be  obtained. 


12.  Muzzle  velocity  expressed  in  terms  of  the  properties  of 
the  propellant 

Assume  that  (m  /M  )«  1  so  that,  in  vievj-  of  Eq.  (32), 

V  V 

Xi  =  . 

Hence  the  muzzle  velocity  per  unit  mass  of  propellant,  namely 
Xi/m^  ,  is  given  by 

&i/m^  = 


or,  in  view  of  Eq,(20),  by 

f 


w 


m 


1  12  }i 


I 


m. 


*'o  «o/ 


RT 


1  ** 


'■^ml 

I  P  I 


-i)a! 


But 


RT/JH  =  (Jr  -D/'d', 
from  Eq.(38),  and  therefore 


m. 


.  /  i 

p,  (2^-1)  A 

/ 

"Ip  / 

(Uo) 


Thus,  the  larger  the  heat  of  combustion  of  the  propellant,  the 
greater  is  the  muzzle  velocity  per  unit  mass  of  propellant. 


13.  Rocket  design  calculations 

The  formulas  developed  in  the  preceding  sections  will  now 
be  applied  to  a  specific  design  problem.  Let  it  be  required 
to  calculate  the  mass  of  the  propellant  and  the  TNT  load  as 
a  function  of  the  chamber  pressure  and  muzzle  velocity  for  a 
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shell  of  diameter  3  in.  and  weight  l5  lb,  exclusive  of  the 
propellant. 

(a)  Calculations  of  muzzle  velocities.  —  The  procedure 
followed  can  best  be  understood  in  connection  with  Table  I. 

'fie  first  assume  that  the  propellant  is  to  be  furnished  with 
an  outside  diameter  of  3A  in.  and  an  inside  diameter  of 
3/16  in.,  so  that  the  web  thickness  D  is  9/32  in.,  or 
0.023it  ft.  A  further  assumption  necessary  is  that  the  maxi¬ 
mum  outside  diameter  of  the  nozzle  is  the  same  as  that  of  the 
shell,  and  that  the  wall  thickness  is  I/I6  in.  The  nozzle 
mouth  is  then  2-7/8  in.  in  diameter  and  the  exit  area  is 
6,ii6  in.^ 

We  assume  a  fixed  mass  m,^  of  propellant,  in  this  case 
3  lb(0.0930  slugs).  The  surface  area  S  as  found  from  3q.(l5a) 
is  369  in.^  Taking  the  chamber  pressure  as  the  independent 
variable  [see  line  (I),  Table  I],  we  find  the  corresponding 
value  of  s/a  from  Fig.  i;.  Since  S  is  369  in.^,  the  corres- 

"D 

ponding  values  of  the  throat  area  A  can  be  found  and  listed 
[line  (3),  Table  I]. 

The  next  step  is  to  set  down  in  line  (U)  the  ratio  of  nozzle 

mouth  area  A  to  throat  area  A. ,  it  being  known  that  A  is 
m  t  m 

6.I46  in.^  Then,  from  Fig.  2,  v/e  find  the  corresponding  pres¬ 
sure  expansion  ratio,  listed  as  (P/Pjjj)^in  line  (5),  Table  I. 

A  comparison  with  line  (2)  shows  that  the  actual  expansion 
ratio  is  appreciably  smaller  than  the  available  pressure  ratio. 
This  is  caused  by  the  fact  that  the  exit  diameter  of  the 
nozzle  is  limited  to  the  shell  diameter.  It  therefore  follows 
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TABLE  I.  Rocket  design  calculations. 


=  3  Ib"*^  S 

=  369  in, 

m 

=  6.U6  in.^ 

D  =0.281 

in. 

(  1) 

P(lb/in.2) 

368 

739 

IU70 

29I4O 

5880 

11,760 

(  2) 

p/p 
'  m 

29 

90 

100 

200 

liOO 

800 

(  3) 

A^(in.2) 

U.61 

2.93 

2.09 

1.62 

1.39 

1.25 

(  h) 

i.Uo 

2.20 

3.17 

3.99 

I4.69 

9.17 

(  5) 

(p/p„) ' 

9.0 

10.  U 

18.2 

26.0 

32.0 

37.0 

(  6) 

Excess  Force 
(lb) 

380 

361 

U27 

636 

1180 

i960 

(  7) 

w^(ft/sec) 

9210 

60it0 

6930 

6820 

6990 

7060 

(  8) 

iTi^(lb/sec) 

11.99 

lii.60 

20.9 

32. 14 

99.6 

99.8 

(  9) 

mdW^(lb)  ! 

1870 

27UO 

iii9o 

6860 

12,000 

21,900 

(10) 

Fsdb) 

2250 

6820 

i;980 

7900 

13,200 

23,900 

(11) 

Wj^(ft/sec) 

6260 

6820 

7210 

7I490 

7690 

7720 

(12) 

Xi;(ft/sec) 

ilUo 

12U0 

1310 

1399 

1390 

lli09 

(13) 

Thrust  Goeff. 

1.32 

l.IiU 

1.92 

1.97 

1.61 

1.63 

(lU) 

\(lb)'" 

2.02 

2.02 

2.86 

U.60 

9.U6 

18. 14 

(19) 

1 — 1 

B 

0.13, 

0.19 

0.26 

0.37 

0.93 

0.7I4 

(16) 

Nozzle  and 
Flange  (lb) 

0.27 

0.93 

0.66 

0.73 

0.76 

0.78 

(17) 

XW(lb) 

!  2.U2 

2.7U 

3.78 

9.70 

10.79 

— 

(18) 

1  (in.) 

11.0 

i 

10.6 

9.31 

6.90 

0.61 

— 

(19) 

TNT  (lb) 

1 

U.72 

14.32 

3.97  . 

1.614 

— 

These  masses  are  expressed  in  pounds  of  mass;  the  values 
may  be  converted  to  slugs  by  dividing  them  by  32.2.  The  mass 
in  pounds  and  the  weight  in  pounds  are  practically  equal 
numerically. 


that  the  nozzle  discharges  the  gases  at  pressures  higher  than 
atmospheric,  and  an  unbalanced  excess  pressure  will  exist 
over  the  nozzle  exit  area.  For  example,  for  P  =  368  Ib/in.^ 
and  a  pressure  expansion  ratio  of  5.0  (see  the  first  column 
of  data  in  Table  I),  the  pressure  is  368/5.0,  or 
73.6  lb/in.2  Taking  the  atmospheric  pressure  as  lit. 7  Ib/in.^, 
the  unbalanced  pressure  is  58.9  Ib/in.^,  and  this,  acting  over 
6.I16  in.^,  produces  an  excess  force  of  38O  lb.  Line  (6), 

Table  I,  gives  the  excess  force  thus  calculated  for  the  other 
pressure  ratios  assumed. 

The  other  component  of  thrust  is  that  caused  by  the 
momentum  of  the  escaping  jet.  The  velocity  can  be  cal¬ 
culated  from  Eq,(20a),  the  value  of  the  radical  being  read 
from  the  curve  in  Fig,  1 .  The  pressure  ratio  to  be  used  is 
that  given  in  line  (5).  The  rate  of  gas  discharge  is 
found  from  Eqs.(7)>  (8)  and  (8a)}  that  is, 

m^  =  /^A^P, 

the  value  of  being  given  in  Sec,  5  as  11,97  slug/lb  sec 
for  the  propellant  used  here.  The  rate  of  discharge  of  the 
gas  is  tabulated  in  line  (8),  and  the  jet  reaction  is  listed 

in  line  (9j. 

The  total  thrpst  Fg  is  obviously  the  sum  of  lines  (6)  and  (9).  From 
F2  the  effective  je'^  velocity  can  be  found}  more  specifi¬ 
cally, 

wl/w„  S'  Column  lO/column  9. 
ra'  m 


o.'? 


z&- 


f/q,  S.  Propa/s/oo  eff/c/eocy  ver>50S  , 
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The  muzzle  velocity  can  now  be  calculated  from  Eq.(31). 
The  initial  mass  of  the  loaded  projectile,  M^,  is  (l3  +3)  lb, 
or  18  lb,  w  being  taken  as  the  effective  velocity  w'  listed 
in  line  (11^  The  muzzle  velocity  is  given  in  line  (12).  As  a 
matter  of  interest,  the  thrust  coefficient  is  given  in  line(13)j 
it  is  calculated  from  Eq.(25a). 

/ie  now  assume  other  values  of  the  propellant  mass  and 
repeat  the  calculations.  The  results  are  shown  in  Figs.  6 
and  8.  Figure  6  is  a  plot  of  the  muzzle  velocity  as  a 
function  of  propellant  mass  (expressed  in  pounds)  for  the 
various  assumed  pressure  ratios.  Figure  8  represents  plots  of 
the  effective  jet  velocity  and  the  thrust  coefficient,  each  as 
a  function  of  pressure  and  two  different  weights  of  propellant 
charge.  Figure  9  gives  the  burning  time  as  a  function  of  the 
pressure  expansion  ratio.  Since  'we  have  assumed  the  same  web 
thickness  in  all  cases,  this  curve  is  the  same  for  all  designs, 
(b)  Calculation  of  weight  of  TMT  charge.  —  The  shell 
structure  is  shown  diagrammatically  in  Fig.  10.  The  shell 
head  is  threaded  on  a  flanged  plate  which  forms  one  end  of 
the  combustion  chamber.  The  cylindrical  portion  of  the  com¬ 
bustion  chamber  is  welded  to  this  plate,  and  the  end  plate 
which  carries  the  nozzle  is  threaded  into  the  chamber  body. 

The  following  data  are  assumed.  The  shell  nose  and  body  are 
made  of  3/l6-in.  steel  and  the  nozzle  is  made  from  steel  of 
thickness  l/l6  in.  The  chamber  and  plates  will  have  thick¬ 
nesses  that  are  determined  by  the  pressure.  The  safe  allovrable 


stress  is  taken  as  $0,000  lb /in.®  in  direct  stress  and 
31, $00  lb/in.2  in  shearing  stress.  The  threaded  flanged  plate 
will  have  a  flange  l/U-in.  long  and  3/l6-in.  thick  as  a  mini¬ 
mum,  unless  the  plate  thickness  is  larger,,  in  which  case  the 
flange  will  be  eliminated.  The  plate  thickness  will  be  calcul¬ 
ated  for  direct  and  shearing  stresses,  and  the  thicker  plate 
will  be  chosen.  The  tvro  plates  will  be  assumed  to  be  of  the 
same  weight. 

The  radius  of  the  shell  ogive  is  taken  as  6  calibers. 

The  surface  of  the  nose  is  UU  in.®,  and  the  volume  is  2$  in,^ 

The  weights  per  unit  volume  of  the  propellant  and  of 

the  TNT  will  both  be  taken  as  0.058  Ib/in.®  The  vreight  per 

unit  volume  of  the  steel  is  taken  as  0.2d  Ib/in.^. 

Shell  nose  and  body.  The  nose  weighs 

UU  X  (3/16)  X  0.28  =  2.31  lb. 

The  shell  body  weighs 

IT  X  3  X  (3/16)  X  0.28  1  =  0.ii9li  1  lb, 

s  s 

where  1  'is  the  length  of  the  shell.  The  volume  available  for 
s 

TNT  is  25  +  J  X  (2  |)  Ig,  or  (25  +  5,. 31  1^)  in.^  The  weight 

of  TNT  is  0^058  (25  +  5.31  1  ),  or  (l.i;5  +  0.308  l  )  Ib.  The 

s  s 

■shell  length  1  will  be  determined  later, 
s 

Weight  of  combustion  chamber.  The  wall  thickness  is  com- 

?/ 

puted  from  Clavarino's  formula,— 

1 

“  ((3^+P)/(3-'r -iiP)]'',  (U) 

7/  For  a  derivation  of  this  formula  see  Martin,  Textbook 
of  mechanics,  Vol.Ili,  ’’Kechanics  of  materials",  pp.  176-180, 
particularly  exercise  2lU, 


-  35  - 


where  Tg  [=  1,5  in.]  is  the  outer  radius  of  the  chamber,  r^  is 
the  inner  radius  and  CT  is  the  stress  (Ib/in. ^).  The  least 
acceptable  thickness  is  assumed  to  be  l/32  in. 

The  length  of  the  chamber,  1  ,  may  be  calculated  from 

c 

the  propellant  weight.  Assume  as  before  that  the  outside  and 
inside  diameters  of  the  propellant  are  3A  in,  and  3/16  in. , 
respectively,  and  that  6  sticks  are  used,  'de  first  check  the 
density  of  loading.  The  net  frontal  area  of  the  6  sticks 
is  2,00  in.^  At  a  pressure  of  800  atmos  the  area  of  the 
combustion  chamber  is  U.50  in.^  Hence  the  density  of  load¬ 
ing  is  2.00/U.50  or  0.1ih5j  which  is  an  acceptable  value.  It 
will  be  somewhat  less  than  this  for  lower  pressures. 


The  length  1  can  then  be  found  from  the  relation 
°  c 

m_  =  6  1  A  , 

r  c  p^ 

where  A^  [=  2.00  in.^]  is  the  cross-sectional  area  of  the  pro¬ 
pellant,  £  [=  0.058  Ib/in.^]  is  the  density  of  the  propellant 


and  m^  is  the  mass  of  the  propellant.  The  chamber  masses  m^ 
as  calculated  from  the  thickness  ~  r^  and  the  length  1^ 
are  given  in  line  (lU^  Table  I. 

End  plates.  To  be  on  the  safe  side,  the  stresses  in  bend¬ 
ing  were  calculated  on  the  assumption  that  the  plates  were 


simply  supported  at  the  edges.  The  maximum  stress  IT  is  at 


the  center  and  is  given  by 

<3-  =1.25  P(r2/tp)%  (U2) 

where  is  the  radius  and  tp  is  the  thickness  of  the  plate. 

For  all  cases  it  turned  out  that  the  bending  stress  was  greater 
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than  the  shearing  stress  and  hence  the  plates  were  designed 
for  bending  rather  than  shear.  The  weights  of  the  two  plates 
are  given  in  line  (l6).  Table  I. 

Nozzle.  The  weight  of  the  nozzle  was  calculated  on  the 
assumption  that  the  wall  thickness  was  l/l6  in.  and  the  cone 
angle  was  l6°.  Allowing  20  percent  increase  in  mass  due  to 
the  flange  and  fillets,  one  may  express  the  weight  in  pounds 
quite  closely  by  the  formula 

=  0.1^  (A^  -  A^).  (U3) 

The  calculated  values  are  given  in  line  (16),  Table  I. 

Weight  of  TNT.  At  this  point  we  can  total  up  all  the 

weights  except  that  for  the  shell  body  and  the  weight  of  TNT 

located  therein.  The  suiji  of  the  chamber,  end  plate  and  nozzle 

weights  is  given  in  line  (17)^  marked  SW. 

To  get  the  total  weighji,  we  add  to  £W  the  weight  of  the 

shell  nose  and  body,  which  is  (2.31  +  0,U9ii  1  )  lb,  and' the 

•  s 

weight  of  TOT,  which  is  l.it^  lb  for  the  TNT  in  the  nose  and 

0,308  1  lb  for  that  in  the  shell  body.  The  grand  total  of 
s 

these  components  must  be  13  lb;  that  is, 

SIW  +  (2.31+0.ii9U  1  )  +  (1.U5+0.308  1  )  =  15  lb, 

0  s 

or 

1^  =  iii.05  -  1.2$  zw. 

The  values  of  1^,  the  leng-th  of  the  shell  body,  are  given 
in  line  (l8j.  Knowing  this,  the  corresponding  total  weights  of 
TNT  can  be  calculated  and  are  given  in  line  (19). 
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Discussion.  —  The  results  of  these  calculations  are  sum¬ 
marized  graphically  in  Fig,  ?•  These  curves  show  the  TNT 
load  as  a  function  of  the  propellant  charge  for  various 
values  of  the  pressure.  The  dotted  lines  are  the  loci  for 
constant  muzzle  velocity.  For  example,  for  every  point  on  the 
dotted  curve  marked  =  1000,  the  corresponding  design  has  a 
muzzle  velocity  of  1000  ft/sec.  These  curves  were  obtained  by 
projecting,  from  Fig.  6,  the  points  of  intersection  of  constant 
X3_-lines  with  constant  P/P^-lines  on  to  the  corresponding  pres¬ 
sure  lines  of  Fig.  ?• 

It  appears  that  the  largest  load  of  TNT  can  be  carried 
when  the  pressure  is  somewhere  between  $0  and  25  atmos.  The 
weight  of  TNT  vdll  then  be  between  U«75  and  5  lb.  For  pres¬ 
sure  either  above  or  below  this  value,  the  "payload"  decreases. 
This  is  to  be  expected.  For  very  low  pressxures,  the  thermal 
efficiency  is  small  and  the  propellant  weight  must  be  increased 
to  attain  a  given  muzzle  velocity.  This  means  a  longer  com¬ 
bustion  chamber.  For  very  high  pressures,  the  thermal  effi¬ 
ciency  increases  very  slowly  but  the  weight  of  the  chamber 
goes  up  rapidly,  and  again  the  payload  will  decrease. 

lU,  Perfonnance  of  geometrically  similar  shells 

These  results  are  more  general  and  therefore  more  widely 
applicable  than  they  may  appear  offhand.  While  the  calcula¬ 
tions  were  made  for  a  shell  of  specific  size,  it  is  easy  to 
show  that  the  conclusions  as  to  optim.um  pressures  apply  to  all 
shells  geometrically  similar  to  one  considered  here. 
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If  the  diameter  and  length  are  multiplied  by  the  same 
scale  factor  s,  the  available  volume  and  v/eight  of  TNT  I’d.!! 
vary  as  s^.  For  constancy  of  energy  (during  fragmentation) 
per  unit  mass  of  shell  nose  or  body,  the  thickness  of  these 
parts  must  therefore  be  increased  by  the  same  scale  factor. 

Hence  the  weight  of  the  shell  and  its  TNT  load  will  vary  as  s^. 

If  all  shells  in  the  series  are  to  have  the  same  muzzle 
velocity  Xi,  it  is  clear  from  Eq,(32)  that  the  mass  of  the 
propellant  charge  must  likewise  vary  as  _s^,  assuming  that  the 
effective  jet  velocity  is  the  same  for  all  shells.  Let  us 
see  what  this  requires.  Assume  tentatively  that  all  shells 
operate  at  the  same  chamber  pressure.  Then,  since  the  throat 
area  A,  and  nozzle  exit  area  A^  both  vary  as  s^,  A  /A^  is 
constant.  Hence  all  nozzles  have  the  sairie  expansion  ratio  and, 
since  the  chamber  pressure  was  assumed  to  be  the  same  for  all 
designs,  it  follows  that  the  effective  jet  velocity  -will  be 
the  same  for  the  entire  series  of  shells.  To  sum  up  to  this 
point,  we  have  the  result  that,  in  a  series  of  geometrically 
similar  shells,  all  working  at  the  same  pressure,  the  muzzle 
velocity  and  jet  velocity  will  remain  constant  provided  that 
the  propellent  weight  is  increased  in  proportion  to  the  weight 
of  the  shell. 

The  argument  is  not  quite  complete  for,  up  to  this  point, 
we  have  ignored  the  weight  variations  of  the  combustion  chamber, 
end  plates  and  nozzles.  From  Eq, (hi),  is  the  same  for 

all  shells  for  equal  pressures  and  equal  allowable  stresses. 


Hence  the  weight  of  the  combustion  chamber,  which  is  propor¬ 


tional  to  l^Crg^-ri^),  must  vary  as  _s^.  Similarly,  an 
examination  of  Eq.(ii2),  which  gives  the  stress  in  the  end 
plates,  shows  that  the  thickness  t^  must  vary  as  s  for 
constant  stress  and  hence  the  weight,  being  proportional 
to  ro^t  must  vary  as  s^.  The  nozzle  weight  will  likewise 

p  — 

be  seen  to  vary  as  s^  if  we  remember  that  Eq.(l4.3)  conceals, 
in  its  numerical  coefficient,  the  length  of  the  nozzle.  This 
length  vail  vary  as  s  provided  that  the  same  taper  is  used  in 
all  designs.  From  the  foregoing  arguments,  it  is  evident  that 
if  the  optimiim  pressure  is  somewhere  betv/een  2$  and  $0  atmos 
for  a  3-in.  shell,  the  same  range  of  pressures  will  be  the 
optimum  for  all  geometrically  similar  shells.  Since  the  optimum 
condition  is  not  critical,  we  can  expect  that  small  deviations 
in  proportions  from  the  ones  assumed  here  will  not  materially 
affect  the  conclusion  as  to  the  best  pressure  range. 

There  is  one  more  point  to  notice.  Constancy  of  pressure 
requires  that  S/A^  be  constant.  Hence  S  varies  as  s^.  But, 
since  the  propellant  mass  m^  varies  as  s^,  the  web  thickness 
D  must  vary  as  s.  Again,  since  S  is  proportional  to  the 
length  of  the  propellant  stick  and  the  mean  diameter  of  the 
stick,  the  latter  must  also  vary  as  s.  From  this  it  follows 
that  the  burning  time  X  will  be  proportional  to  the  scale 
factor  but  that  the  niomber  of  sticks  of  propellant  is  to  be 


kept  the  same  for  all  designs. 
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APPENDIX 

List  of  Symbols  Used 

Area  of  nozzle  mouth  (ft®) 

Cross-sectional  area  of  propellant  (in.^) 

Area  of  nozzle  throat  (ft^) 

Propellant  burning  constant  (ft/sec) 

Propellant  burning  constant  [ (ft/sec)/(lb/ft^) ] 
Design  parameter 

Mean  specific  heat  at  constant  pressure  (Btu/slug  °F) 
Mean  specific  heat  at  constant  volume  (Btu/slug  °F) 
Web  thickness  of  propellant  (ft) 

Efficiency 

Thrust  for  simple  nozzle  (lb) 

Thrust  for  expanding  nozzle  (lb) 

Acceleration  due  to  gravity  (ft/sec®) 

Heat  of  combustion  of  propellant  (Btu/slug) 

Mechancial  equivalent  of  heat  [=  778  ft  Ib/Btu] 

Mass  of  gas  discharged  per  unit  time,  per  unit  of 
pressure  [ (slug/sec)/(lb/ft2) ] 

Length  of  chamber  (ft) 

Length  of  shell  body  (ft) 

Mass  of  propellant  charge  (slug) 

Mass  burned  per  unit  time  (slug/sec) 

Mass  of  gas  in  chamber  (slugs) 

Mass  of  gas  discharged  per  unit  time  (slug/sec) 


Instantaneous  mass  of  projectile  (slug) 

Initial  nraSIs  of  projectile,  including  propellant 

ch^ge  (slug) 

Initial  mass  of  projectile,  exclusive  of  pro¬ 
pellant  charge  (slugs) 

Chamber  pressure  (Ib/ft^) 

Atmospheric  pressure  (Ib/ft^) 

Bomb  pressure  (Ib/ft^) 

Equilibrium  chamber  pressure  (Ib/ft^) 

Pressure  at  nozzle  mouth  (Ib/ft^) 

Pressure  at  nozzle  throat  (Ib/ft^) 

Gas  constant  for  an  ideal  gas  (ft  Ib/slug  °F) 

Inner  radius  of  combustion  chamber  (in. ) 

Outer  radius  of  combustion  chamber  (in.) 

Total  burning  surface  (ft  ) 

Time  (sec) 

Thickness  of  end  plates  (in. ) 

Temperature  of  gas  (°F  absolute) 

Specific  volume  of  gas  (ftVsl^g) 

Instantaneous  volume  of  gas  in  chamber  (ft^) 
Initial,  or  clearance, volume  (ft^) 

Initial  voliime  of  propellant  (ft^) 

Gas  exit  velocity  (ft/sec) 

Effective  gas  exit  velocity  (ft/sec) 

Gas  velocity  at  nozzle  throat  (ft/sec) 
vYeight  of  nozzle  and  flange  (lb) 

Velocity  of  projectile  at  time  t 
Muzzle  velocity  of  projectile  (ft/sec) 
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Z 

i 

X 

Xe 

Ac 

T 


Fractional  mass  of  propellant  burned 
Propellant  constant 
M  /m 

0'  T 

Propellant  constant  (slug/lb  sec) 

Adiabatic  exponent  for  propellant  gases 
Density  of  propellant  (slug/ft^,  Ib/in.^) 
Ratio  of  chamber  pressure  to  bomb  pressure 
Equilibrium  pressure  ratio 
Initial  pressure  ratio 
Stress  (Ib/in.^) 

Burning  time  (sec) 


